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neighbor  0-0  distance  ( roo^)  between  al)out  2.5  and  2.9  A.  However,  the  data 

can  also  be  fitted  with  a  slight  cubic  dependence,  which  when  extrapolated, 
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gives  rise  to  a  minimum  in  Av  at  i*QQ  =  2.35  ±  0.06  A.  This  is  close  to  the 

A 

value  of  2.4  A  found  for  the  symmetric  O-H-O  units  in  other  material s^This 
minimum  suggests,  but  does  not  prove,  that  a  symmetric  hydrogen  bond  struc¬ 
ture,  O-H-O,  and  thus  a  new  ice,  may  result  at  pressures  of  about  75±20  GPa. 
Treatment  of  the  0-H***0  units  according  to  coupled  oscillators  using  experi- 
mentally  confirmed  assumptions  indicates  that  Av  ^  A/rQH,  which  when 

combined  with  present  and  previously  published  data  yields,  approximately, 

•v  -1°2 

A  >  2943  cm  A  .  Thus  the  decrease  observed  in  Av  is  thought  to  occur 
because  the  OH  bond  length  ( r0H^  increases  increasing  pressure  as 

the  0-0  distance  rnn  in  the  corresponding  hydrogen-bonded  0-H***0  unit 


decreases. 


RAMAN  AND  X-RAY  INVESTIGATIONS  OF  ICE  VII  TO  36.0  GPa 
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and 

F.  A.  Mauer,  S.  Block,  G.  J.  Piermarini,  and  R.  Munro 
Center  for  Materials  Science 
National  8ureau  of  Standards 
Washington,  DC  20234 

I.  INTRODUCTION 

Ice  VII,  discovered  by  Bridgman'*’,  is  produced  when  ice  VI  is  subjected 

2 

to  a  pressure  £  2.14  GPa  at  room  temperature  .  X-ray  diffraction  data  at  high 

pressure  indicate  that  ice  VII  has  a  body-centered  cubic  arrangement  of  oxygen 

3  4  5 

atoms,  with  some  degree  of  disorder  of  the  protons  ’  ’  .  In  contrast,  ice 

VIII,  which  forms  when  ice  VII  is  cooled  below  -3  °C,  is  proton-ordered6, 

and  its  oxygen  arrangement,  although  originally  thought  to  be  the  same  as 

that  of  ice  VII,  probably  involves  slight  tetragonal  distortion,  space 
7  8 

group  14^/amd  ’  . 

An  interesting  feature  common  to  ices  VI,  VII,  and  VIII,  involves 
interpenetration  of  lattices.  In  the  case  of  ice  VII  two  interpenetrated, 
but  unconnected,  diamond-type  (like  ice  Ic)  lattices  exist.  Thus,  the 
density  of  ice  VII  approaches  twice  that  of  ice  Ih,  and  would  be  twice  if 

It  were  not  for  an  increase  in  the  0-0  distance  rgg  above  the  ice  Ih 

e 

value  of  2.76  A.  A  further  feature  of  ice  VII  that  is  related  to  the 

unconnected  lattices,  is  that  only  four  of  the  eight  nearest-neighbor 

oxygen  atoms  are  hydrogen-bonded  to  the  central  oxygen  atom.  The 

oxygen  atoms  of  the  two  unconnected  lattices  engage  in  repulsive 

o 

contacts  which  increase  rgg  above  2.76  A. 
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Of  all  of  the  phases  of  water  known  to  exist  thus  far,  the  fluid  phase 
(including  the  liquid  and  the  dense  supercritical  fluid),  and  ice  VII,  have 

g 

the  greatest  pressure  range  of  stability  .  Ice  VII  has  been  visually  observed 

with  a  polarizing  microscope  to  50  GPa  at  room  temperature  with  no  obvious 

g 

transition  and  recently  confirmed  to  36  GPa  by  x-ray  powder  diffraction 
measurements1®.  No  ice  having  symmetric  hydrogen  bonds  has  yet  been  produced. 
Such  an  ice  would  require  very  high  pressures  for  its  formation,  calculated 
to  be  between  35  and  80  GPa11. 

The  extreme  pressure  stability  of  ice  VII  is  a  feature  which  made  the 

present  Raman  and  x-ray  study  attractive.  This  large  pressure  range  allowed 

for  large  variations  in  Av  and  in  the  x-ray  lattice  parameter,  a,  compared, 

12 

for  example,  to  ice  VI,  where  a  similar  study  was  impracticable  .  Between 
2.30  GPa  and  36  GPa  the  body-centered  cubic  unit  cell  edge,  a,  decreased 

O 

from  3.367  to  2.896  A,  and  to  30  GPa,  the  Raman  OH  peak  frequency,  Av, 
decreased  from  about  3300  to  2800  cm'1.  These  large  decreases  allowed  the 
pressure  dependences  to  be  determined  accurately,  which  in  turn  allowed  for 
detailed  analyses  and  interpretations  of  the  data. 

II.  EXPERIMENTAL 
X-ray  Method 

Measurements  of  the  pressure  dependence  of  the  lattice  parameter  of 
Ice  VII  were  made  by  the  energy  dispersive  powder  diffraction  method  using 
a  diamond  anvil  pressure  cell  .  The  ruby  fluorescence  method  was  used  to 
measure  the  pressure1^ ,1^.  The  pressure  cell  was  mounted  on  a  horizontal 
diffractometer  equipped  with  a  tungsten  tube,  an  intrinsic  germanium  detec¬ 
tor,  and  a  0.15  mm  i.d.  collimator  made  from  telescoping  sections  of  tungsten 
carbide  tubing.  An  x,  y,  z  micrometer  stage  was  used  for  mounting  and  - 
aligning  the  cell  with  respect  to  the  x-ray  beam  and  the  diffractometer 


axis.  With  proper  alignment,  patterns  could  be  recorded  with  no  trace  of 

diffraction  lines  from  the  gasket.  The  sample  was  also  centered  on  the 

axis  of  rotation  so  that  it  could  be  oscillated  through  ±  1.7°  to  reduce  the 

effect  of  sample  graininess,  a  severe  problem  because  of  the  small  sample 
-3  3 

volume  (2  x  10  mm  ). 

In  addition  to  the  0.2°  receiving  slit,  a  0.15  mm  wide  slit  on  the 
detector  arm,  2.5  cm  from  the  diffractometer  axis,  was  used  to  exclude  from 
the  detector  the  radiation  scattered  along  the  path  of  the  primary  beam 
through  the  diamonds.  This  slit  reduced  the  background  counts  by  about  a 
factor  of  ten  over  much  of  the  energy  range  and  greatly  improved  the  signal- 
to-noise  ratio. 

A  diffraction  angle  of  20  =  14.00°  was  selected.  Assuming  that  the 
usable  range  of  the  x-ray  continuum  extends  from  10  keV  (below  which  less 
than  2  percent  of  the  incident  flux  is  transmitted  by  the  diamonds)  to 

o  o 

50  keV,  the  range  of  possible  d  spacings  extends  from  5.1  A  down  to  1.0  A, 
equivalent  to  the  range  obtained  with  a  20  scan  from  17°  to  99°  with  CuKa 
radiation. 

Generally  the  lattice  parameter  at  each  pressure  was  based  on  the  (110), 
(200),  and  (211)  reflections.  The  precision  of  the  lattice  parameter  measure¬ 
ment  was  about  0.25  percent  while  the  accuracy  was  about  1  percent.  A  more 
significant  error  resulted  from  the  uncertainty  of  the  pressure  in  the  highest 
ranges  arising  from  the  nonhydrostatic  character  of  ice  VII.  For  example, 
the  pressure  measured  by  the  ruby  fluorescence  method  in  a  sample  approximately 
0.2  mm  in  diameter  by  0.1  mm  thick  ranged  from  23.7  GPa  to  26.2  GPa.  This  range 
corresponds  to  0.6  percent  uncertainty  in  a  and  0.8  percent  uncertainty  in  Av. 

To  circumvent  this  difficulty,  the  average  pressure  measured  by  ruby  chips 
situated  throughout  the  sample  was  employed. 
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Raman  Method 

Raman  spectra  of  ice  VII  were  obtained  using  the  diamond  anvil  pressure 

cell  and  the  ruby  fluorescence  method  for  pressure  measurement.  The  spectra 

were  obtained  using  two  different  double  monochromators.  For  pressures  to 

* 

15  GPa  a  standard  Instruments  S.A.  HG-2S  instrument  was  employed,  and  for 
pressures  from  15  to  30  GPa  an  HG-2S  instrument  equipped  with  red-sensitive 
holographic  gratings  was  used.  Detection  was  accomplished  with  cooled  and 
uncooled  Hamamatsu  R928  (S-20)  PM  tubes,  a  Fluke  412B  HV  supply,  a  Keith! ey 
414S  pi  coammeter,  an  Esterl ine- Angus  L1101S  high-speed  recorder,  and  at  the 
highest  pressures  with  a  Nicolet  multichannel  analyzer.  Slit-widths  cor¬ 
responding  to  16  cm  ^  were  used  in  all  cases,  with  scanning  rates  of 
50  cm  Vin  1.  Excitation  was  accomplished  with  514.5  and  488.0  nm  excitation 
from  Spectra  Physics  model  165  argon  ion  lasers.  At  pressures  from  15  to 
30  GPa,  an  interference  filter  was  employed  with  488.0  nm  excitation. 

The  diamond  anvil  cell  was  used  in  two  orientations.  To  pressures  of 
15  GPa,  scattering  at  an  acute  angle  was  employed  with  collection  in  the 
forward  direction^-  .  At  pressures  from  15  to  30  GPa  the  acute  angle  scattering 
was  also  used  but  collection  was  accomplished  in  the  backward  direction. 
Pressures  were  measured  simultaneously  with  Raman  spectra  by  making  compari¬ 
sons  between  the  fluorescence  from  the  ruby  in  the  cell,  and  from  an  external 
ruby  in  each  run.  Light  scattered  from  an  unfocussed  laser  beam  at  a  power 
level  of  100  mW  was  used  in  pressure  measurements  and  care  was  taken  not  to 
focus  on  the  ruby  chips.  In  contrast  a  sharp  focus  was  used  in  the  ice  to 
reduce  the  scattering  from  the  diamonds  during  Raman  measurements.  No 

^Certain  commercial  equipment,  instruments,  or  materials  are  identified  in 
this  paper  in  order  to  adequately  specify  the  experimental  procedure.  Such 
identification  does  not  imply  recommendation  or  endorsement  by  the  National 
Bureau  of  Standards,  nor  does  it  imply  that  the  materials  or  equioment  iden¬ 
tified  are  necessarily  the  best  available  for  the  purpose. 
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marked  broadening  and  loss  of  detail  in  the  two-phonon  diamond  band  was 
observed  in  this  work  to  30  GPa. 

For  pressures  up  to  20  GPa,  the  precision  of  the  Raman  measurements  was 
about  2  cm”*  while  the  accuracy  was  5  cm”*.  Above  20  GPa,  the  accuracy  was 
10  cm”*. 

III.  EXPERIMENTAL  RESULTS 
X-ray  Oata 

The  lattice  parameters  obtained  for  ice  VII  to  36  GPa  are  tabulated  in 

Table  I,  and  shown  versus  pressure  in  Fig.  1.  From  Fig.  1  it  is  apparent 

o 

that  the  lattice  parameter,  a,  (in  A)  decreases  at  a  decreasing  rate  with 
pressure  rise.  Various  polynomials  giving  pressure  P  as  a  function  of 
lattice  parameter  were  fit  to  these  data.  Of  these,  a  cubic  fit  using  a 
digital  computer  with  72-bit  words  reproduced  the  data  within  4  percent. 

The  equation  with  pressure  in  GPa  and  lattice  parameter  in  angstroms  was 
as  follows: 

P  =  7.8166  x  103  -  6.8985  x  103  (a)  +  2.0368  x  103  (a)2  -  2.0114  x  102(a)3.  (la) 

(The  coefficients  are  given  to  five  figures  to  represent  the  cubic  equation 
accurately.  These  numbers  should  not  be  considered  independently  as  signifi¬ 
cant  figures.)  The  least  squares  fit  is  shown  by  the  solid  line  through  the 
data.  Fig.  1. 

Volumes  determined  from  the  present  measurements  differ  by  about  2.3  per¬ 
cent  from  Bridgman's  early  work*^.  There  is  good  agreement  with  the  results 

18  19 

of  Olinger  and  Halleck  ,  when  their  revised  pressure  calibration  is  used  . 

20 

An  absolute  comparison  with  the  results  reported  by  Holzapfel  and  Drickamer 
could  not  be  made  because  their  absolute  volume  measurements  were  not  given, 
but  a  relative  comparison  of  V/VQ  shows  good  agreement  out  to  about  12.5  GPa. 


6 


The  divergence  observed  at  higher  pressures  can  be  attributed  largely  to  the 
different  pressure  calibrations  used  which  differ  by  about  2.6  percent  at 
the  higher  pressures. 

Raman' Data 

Raman  OH-stretching  peak  frequencies  obtained  for  ice  VII  to  30  GPa 
are  tabulated  in  Table  II,  and  plotted  in  Fig.  2.  The  peak  frequency,  Av, 

(in  cm  *)  decreases  at  a  decreasing  rate  with  pressure  rise.  This  decrease 
is  similar  to  that  observed  for  the  x-ray  data. 

The  Raman  data  were  again  treated  by  least  squares  polynomials.  The 
fit  of  a  cubic  equation  which  reproduced  the  data  within  3.8  percent  with  P 
as  a  function  of  Av,  was  as  follows: 

P  =  6.7950  x  103  -  6.2337  (aG)  +  1.9164  x  10"3  (Au)2  -  1.9732  x  10'7  (aG)3.  (lb) 

(Again  P  is  in  GPa,  aG  is  in  cm  1,  and  the  five  digits  do  not  independently 
represent  significant  figures.)  The  least  squares  fit  is  shown  by  the  solid 
line.  Fig.  2. 

Raman  spectra  typical  of  those  obtained  in  this  work  in  the  OH-stretching 
region  are  shown  in  Fig.  3  for  pressures  of  3.2  GPa  and  15.5  GPa.  In  addi¬ 
tion  to  the  obvious  decrease  in  the  frequency  of  the  Raman  OH-stretching 
peak  position,  the  shoulders,  evident  at  3.2  GPa  between  about  3350  cm  1  and 
3600  cm  *,  also  show  frequency  decreases  with  pressure  rise.  However,  the 
shoulder  intensities  wane  relative  to  the  peak,  making  detection  very  diffi¬ 
cult  above  about  10  to  15  GPa.  When  the  OH-stretching  peak  is  distinct  from 
the  shoulders,  e.g. ,  at  15.5  GPa,  Fig.  3,  the  half-width  (full  width  at  half¬ 
height)  of  the  peak  corrected  for  slit  width  may  be  estimated,  viz.  ,  85  ± 

15  cm"1. 


-5W  -. 
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Correlation  of  Raman  and  X-ray  Results 
Method  A- Li  near  Fit 

The  x-ray  equation  giving  P  as  a  function  of  a  and  the  Raman  equation 

giving  P  as  a  function  of  Av  were  solved  using  Tqq  =  aV^2  to  yield  values  of 

Av  versus  r^g  for  uniformly  selected  values  of  Tqq  in  the  observed  pressure 

range.  To  examine  the  functional  dependence  of  Av  on  rgg,  polynomials  in 

r go  were  fit  to  the  derived  set  of  points.  A  first  degree  polynomial, 

-  2  3 

Av  =  -5.0725  x  10  +  1.3135  x  10  Tqq,  reproduced  the  data  points  to  0.2  percent. 

-1°-1 

The  slope  of  this  linear  relation,  1313.5  cm  A  ,  is  not  too  far  from  the 

-1°1  21 
value  1520  cm  A  reported  for  ice  1^  at  255  °K  on  the  basis  of  Raman  data  . 

In  Fin.  4,  the  derived  data  points  and  the  straight  line  fit  (curve  A)  are 

shown.  The  gross  trend  of  the  data  is  clearly  represented  by  the  straight 

line  in  the  observed  range.  However,  extrapolation  of  the  linear  fit  to 

o 

small  rQ0  values,  e.g. ,  to  2.4  A  or  less,  presents  a  serious  difficulty. 

Such  a  linear  extrapolation  to  rgg  values  below  the  value  corresponding 
to  a  symmetric  hydrogen  bond,  0-H-0,  is  thought  to  be  equivalent  to  the 
physically  erroneous  situation  in  which  the  OH  stretching  force  constant 
decreases  as  the  OH  bond  strengthens  by  virtue  of  bond  shortening.  Hence, 
a  cubic  fit  was  tried,  because  the  cubic  extrapolation  avoids  a  grossly 
false  physical  situation. 

Method  B— Cubic  Fit 

o 

Data  points  for  Tqq  values  above  about  2.7  A  were  observed  to  deviate 
systematically  from  a  straight  line.  This  region  involves  the  highest  accu¬ 
racy  of  the  Raman  and  x-ray  data,  as  well  as  pressure  measurements,  and  it 
refers  to  a  region  of  complete  hydrostaticity.  No  method  of  data  presenta¬ 
tion,  including  plots  of  raw  data,  has  thus  far  removed  this  curvature.  -To 
further  test  the  validity  of  this  nonlinearity,  the  deviation  from  a  straight 


line,  compared  to  the  uncertainty  in  Au  viewed  as  a  function  of  Tqq,  was 
examined.  An  estimate  of  the  uncertainty  at  the  95  percent  confidence  limit 
yields  uncertainties  of  0.056  percent  for  0  <  P  <  10  GPa;  0.98  percent  for 
10  <  P  <  17.5  GPa;  and  0.20  percent  for  17.5  GPa  <  P  <  30  GPa.  On  this 

basis,  the  straight  line  is  an  adequate  representation  of  the  details  only 

for  the  highest  pressure  range,  i.e.,  the  smallest  values  of  ^g.  The 

largest  excursion  from  a  straight  line  occurs  in  the  lowest  pressure  range 

o 

at  about  rgg  =  2.785A  with  a  deviation  of  about  0.24  percent  which  is  greater 
than  the  inherent  uncertainty  by  more  than  a  factor  of  4.  Therefore,  for 

detailed  accuracy,  a  higher  order  polynomial  should  be  considered. 

A  quadratic  equation  reproduces  the  data  points  only  to  0.15  percent, 
which  is  still  not  adequate.  It  also  gives  rise  to  a  maximum  at  low 
pressure,  and  hence  cannot  also  yield  a  minimum  at  high  pressure.  Therefore, 
its  high  pressure  extrapolation  involves  the  same  problem  mentioned  pre¬ 
viously  with  regard  to  the  linear  extrapolation.  A  cubic  equation,  however, 
successfully  describes  the  data  to  0.03  percent  which  is  compatible  with  the 
uncertainty  of  Au  as  shown  by  curve  B  in  Fig.  4.  The  cubic  equation  is 

Au  =  7.8053  x  104  -  8.6353  x  104  rQQ  +  3.2573  x  104  rQQ2  -  4.0297  x  103  r003‘ 
The  cubic  fit  also  has  the  advantage  that  it  yields  a  realistic  high 

pressure  extrapolation.  At  very  high  pressures  a  symmetric  hydrogen  bond, 

0-H-0,  is  expected  to  occur,  and  this  would  produce  a  minimum  in  the  Au  vs 

Tqq  curve.  The  cubic  equation  thus  allows  an  estimate  of  the  symmetric 

hydrogen  bond  length  to  be  made  by  finding  the  minimum  which  occurs  in  the 

extrapolation  to  smaller  rgg.  The  cubic  equation  yields  rgg(min)  = 

O  - 

2.35  ±  0.06A,  using  a  very  conservative  high  pressure  uncertainty  in  Au  of 
±  20  cm  *.  Also  using  the  pressure  relation,  P  vs  a,  a  corresponding  value 
for  P  and  its  uncertainty  can  be  determined.  P[rgg(min)]  =  75±20  GPa. 
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This  range  of  pressure  agrees  with  a  previous  estimate  of  35  to  80  GPa^ 

O 

and  the  ^(min)  value  agrees  with  a  value  of  2.4  A  presented  subsequently. 
Mode  Gruneisen  Parameter 

The  mode  Gruneisen  parameter  is  a  nonthermodynamic  quantity  that  can 
be  determined  from  the  volume  dependence  of  a  given  vibrational  (Raman  or 
Infrared)  frequency.  The  definition  for  the  ith  Raman  mode  Av.  is  as 
follows: 


3£nAv . 

Y1  32nv  *  ^ 

where  v  is  the  cell  volume,  may  be  either  positive  or  negative,  but  the 
total  thermodynamic  quantity,  y,  where 

Y  =  2  Yi  ,  (3) 

i 

is  positive  for  most  systems  as  seen  from, 

Y  =  v(3p/3T)v/Cv  =  v(3S/3V)t/Cv.  (4) 

Values  of  YqH  were  calculated  for  the  range  of  r^g  values  from  2.35  to 

o 

3.035  A,  using  both  the  method  A  and  method  B  fits  for  Av  versus  Tqq.  The  results 
are  shown  in  Fig.  5. 

Values  from  Method  A  Fit 

A  plot  of  In  Av  vs  In  a  was  found  to  be  roughly  linear,  and  to  have  a 
positive  slope.  Hence,  Yq^  *  constant,  from  Eq.  (2).  To  determine  Yqh  as  a 
function  of  rgg,  the  linear  fit,  method  A,  Fig.  4,  was  employed.  This  pro¬ 
cedure  yielded  the  result  labeled  A  in  Fig.  5.  The  Yq^  ^rom  method  A  is 
roughly  constant,  but  it  shows  a  slight  nearly-1 inear  decrease  with  decreasing 

O 

Tqq.  In  the  region  of  data  between  2.5  and  2.9  A,  the  average  value  of  YqH 
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Is  -0.39.  This  is  close  to  the  Yg^  of  "0.36  reported  for  the  OH  stretching 
23 

■ode  of  ice  1^  .  However,  the  decreasing  trend  of  Yg^  with  decreasing  Tqq 

leads  to  the  same  unrealistic  physical  situation  mentioned  previously  for 
the  linear  extrapolation  of  Fig.  4. 

Yq^  Values  from  Method  B  Fit 

In  Fig.  5,  Yq^  values  determined  from  the  method  B  cubic  fit  of  the 
Fig.  4  data  are  shown,  labeled  B.  Here  a  quadratic  dependence  between 

O 

Yq^  and  Tqq  is  evident  which  yields  a  minimum  near  2.7  A.  The  average  value 

o 

of  Yqh  in  the  region  of  data  between  2.5  and  2.9  A  is  -0.36  in  agreement 

23 

with  the  value  of  -0.36  reported  for  ice  Ih  .  The  extrapolation  of  YqH 

O  O 

to  large,  3.035  A,  and  small,  2.35  A,  values  of  rgg  using  method  B  is  shown 
by  the  dashed  line  in  Fig.  5  (also  for  method  A).  The  high  pressure 
extrapolation  suggests  that  Yg^  rosy  become  positive  at  pressures  where 
symmetric  O-H-O  units  occur.  The  reasons  for  this  are  described  subsequently. 

The  negative  values  of  yqh  obtained  either  by  method  A  or  B  imply  from 
Eq.  (4)  that  other  vibrational  modes  must  have  positive  y^  values.  Ice  1^, 

23 

for  example,  has  a  y.,-  value  of  1.86  for  the  intermolecular  0-0  stretching  mode  , 
whereas  the  intramolecular  OH-stretching  vibration  has  a  YqH  of  -0.36.  A  similar 
situation  is  expected  to  occur  for  ice  VII  at  low  pressures. 

IV.  INTERPRETATION 

Relationship  Between  TqH  and  r^g 

An  approximate  treatment  leading  to  an  inverse  square  relationship 

between  ^  and  Ao  follows.  This  approach  treats  a  hydrogen  bond  in  ice 

VII  as  a  weakly  coupled  oscillator  system,  and  refers  specifically  to  the 

o 

weak-coupling  region  2.54  <  rgg  <  2.90  A.  The  treatment  does  not  apply  to 
values  of  rgg  «  2.54,  where  the  coupling  between  0-H  and  H****0  of  the 
0-H****0  bonds  probably  becomes  extremely  strong. 


o 


Consider  a  hydrogen  bond,  0-H****0,  as  a  weakly  coupled  oscillator 


system  having  a  force  constant  k^  for  the  OH  bond  and  a  force  constant  k^ 
for  the  weak  liaison.  The  system  is  pictured  schematically  in 

Fig.  6.  In  the  0-H****0  linkage  the  OH  distance  is  rg^,  the  H****0  distance, 
r^,  and  rgg  =  rgH  +  r^g.  If  the  magnitude  of  the  proton  displacement  from 
Its  equilibrium  position  is  d,  the  absolute  magnitude  of  the  harmonic 
restoring  force  opposing  this  displacement  is  given  by, 

F  =  kxd  +  k2d,  (5) 

or  by 


l^d  X2d 

r0H  rHQ 


(6) 


24 

where  and  X2  are  the  elastic  moduli  corresponding  to  the  OH  and  H---*0 
bonds,  respectively.  For  convenience,  let  F’  =  F/d,  whence, 


I 


F' 


(7) 


Next,  it  is  assumed  that  X^  »  X2,  l.e. ,  weak  coupling.  Further  it  is 
assumed  that  X^  is  constant,  i.e.,  t  f(rgg).  (This  assumption  may  not  be 
exact,  but  it  has  the  effect  of  forcing  all  of  the  change  in  to  occur 

o 

via  a  change  in  rg^.  Also,  the  increase  in  rgH  for  2.54  <  Tqq  <  2.90  A  is 
estimated  to  be  only  about  9  percent.)  Equation  6  is  thus  reduced  to 


F'  - 


and 


(8) 
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differentiation  with  respect  to  r^  gives, 

dF'  _  .  *1  dr0H 

- - T—  -  (9) 

dr00  r  OH  dr00 

2 

Further  for  a  harmonic  oscillator,  F'  *  cur,  where  c  is  a  constant,  and  to 
Is  the  vibrational  frequency  in  cm"*.  Differentiation  gives  dF'  =  2aodtu, 
and  substitution  in  Eq.  (8)  yields, 


2 audio  „  *1  dr0H 

-  ~  -  — - .  (1 

dr00  r  OH  dr00 

Taking  cfco/drQ0  approximately  constant  =  0,  as  seen  from  A  in  Fig.  4,  and 
assuming  also  that  drQ>/drQ0  =  ~a  -  constant,  then 


or, 


Av 


(ID 

(12) 


where  u>  =  Av,  and  A  =  aK^/2cp. 

Equation  (12)  has  been  tested  in  the  following  manner  assuming  that 
the  relationship  between  rgg  and  r^  is  the  same  whether  deduced  for 
(A)  systems  of  differing  chemical  elements  at  one  pressure,  or  (B)  a  given 
system  in  which  changes  in  bond  distances  are  induced  by  changes  in  pressure. 
In  recent  work,  a  plot  of  rg^  versus  rgg  for  various  hydrogen-bonded 
systems,  at  one  atmosphere  of  pressure  and  having  O-H-O  angles  greater  than 
160®,  was  given  using  carefully  selected  neutron  and  x-ray  data.  Data 
from  that  work  have  been  reproduced  in  Fig.  7.  In  the  region  of  interest, 

O 

2.54  <  Tqq  <  2.90  A,  the  approximation  of  dr^/dr^g  =  contant  =  -a  is  satis¬ 
fied  with  a  =  0.24.  For  each  rgg  in  Fig.  7,  a  value  of  Av  was  determined 
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from  the  cubic  equation  used  In  Fig.  4  which  was  determined  for  one 

system  (H^O,  Ice  VII)  under  pressure.  This  gave  a  new  set  of  values,  Av 

02 

versus  r^,  to  which  Eq.  (12)  was  fit.  The  value  A  =  2943  cm  A  was 

obtained,  and  this  produced  the  solid  curve  shown  In  Fig.  7.  Additionally, 

It  can  also  be  seen  from  Fig.  7,  as  obtained  from  Ref.  (25),  that  the 

o 

symmetric  O-H-O  unit,  i.e. ,  r^  =  2^  occurs  at  2.4  A  for  other  materials. 

O 

This  compares  very  favorably  with  our  own  estimate  of  2.35  ±  0.06  A. 

Equation  (12)  and  its  derivation,  albeit  approximate,  lead  to  implica¬ 
tions  of  relevance  to  all  the  results  of  this  work. 


When  pressure  is  applied  to  ice  VII,  a,  and  thus  rgg,  decrease  at  a 
decreasing  rate.  Fig.  1.  Pressures  in  excess  of  36  GPa  would  only  be  expected 
to  produce  a  further  decrease  in  rgg,  with  the  rate  of  decrease  continuing  to 
decline.  No  minimum  in  the  curve  of  rgg  versus  pressure  would  be  expected. 

The  plot  of  Av  versus  pressure,  Fig.  2,  has  a  curvature  similar  to 
that  of  Fig.  1.  In  contrast  to  the  x-ray  data,  however,  a  minimum  in  Av 
would  be  expected  to  occur  at  some  very  high  pressure,  e.g.  ,  as  shown  by 
curve  B  in  Fig.  4.  Further,  Fig.  2  is  thought  to  indicate  that  the  0H- 
stretching  force  constant  decreases  with  pressure  rise  because  the  OH  bond 


Is  lengthening.  However,  a  pressure,  P  ,  will  eventually  be  obtained  at 
which  =  Tqq/2.  Above  this  pressure,  i.e.,  when  P  >  P$,  the  OH  distance 
may  decrease  at  the  same  rate  as  the  0-0  distance  decreases  (provided  that 


rg^  continues  to  equal  rgg/2) .  If  the  OH  bond  distance  decreases,  the  en¬ 
force  constant  will  increase,  and  Av  will  increase.  Hence,  a  minimum  in  Av 


versus  P,  or  in  Av  versus  rgg,  is  expected,  Fig.  4  curve  B. 

Associated  with  the  occurrence  of  a  minimum  in  Av  vs  rgg  would  be  a 
change  of  sign  for  the  mode  Gruneison  parameter  YqH*  This  would  occur  because 
for  P  >  Pg,  Av  would  increase  with  Increasing  pressure  while  the  volume 
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would  continue  to  decrease.  Hence,  by  Eq.  (2),  Yqh  would  be  positive  for 

P  >  Pg.  At  room  pressure,  Yq^  is  negative  and  decreases  with  an  increase  of 

pressure.  Consequently,  the  hypothesis  of  a  symmetric  hydrogen  bond  at 

high  pressure  and  the  values  of  YgH  at  low  pressure  imply  that  YqH  vs  rQQ 

has  a  minimum  negative  value  and  rises  to  zero  at  the  high  pressure  symmetric 

bond  length.  The  values  of  YqH  computed  by  method  B  and  their  smooth  curve 

extrapolation,  which  are  shown  in  Fig.  5,  exhibit  this  behavior. 

The  minimum  of  YqH  vs  rgg  may  possibly  be  understood  in  terms  of  the 

approaching  dominance  of  the  repulsive  part  of  the  potential.  It  is  known, 

for  example,  that  y  generally  increases  with  density  when  the  repulsive  part 

26 

of  the  potential  behaves  like  that  for  hard  spheres  .  This  situation  may 

be  analogous  to  that  of  symmetric  ice  VII.  This  ice  may  no  longer  contain 

discrete  H2O  molecules  and  may  be  ionic  and  proton  ordered  as  the  result  of 

symmetric  O-H-O  units.  This  situation  is  possible  because  the  r^  distance 

increases  with  pressure  rise,  and  thus  the  s-p  orbital  overlap  between  H  and 

0  may  decrease,  leading  to  higher  ionic  OH  character  compared  to  isolated 

H2O  molecules.  Also,  extremely  strong  0-0  repulsive  forces,  resembling  those 

between  hard  spheres  at  high  density,  would  be  present. 

Finally,  other  features  of  the  present  results  are  of  some  interest. 

-  -1  0 

Figure  4  shows  a  minimum  in  Av  at  2712  cm  for  Tqq  =  2.35  A,  and  a  maximum 
-  -1  0 

in  Av  of  3355  cm  at  rgg  =  3.04  A.  The  pressure  corresponding  to  the 
minimum  Av  is  estimated  to  be  75  ±  20  GPa.  Extrapolation  of  the  cubic  relation 
to  the  maximum  is  not  meaningful  because  the  maximum  occurs  in  a  physically 
unobservable  region  .  In  regard  to  this,  the  Raman  data  of  Fig.  2  below  5  GPa 
were  treated  by  a  first  degree  polynomial  with  Av  as  a  function  of  P.  For 
P  *  0,  which  essentially  corresponds  to  the  maximum  feature  in  Fig.  4,  Av  = 

3378  cm  *.  The  value  of  3378  cm  *  may  refer  to  an  upper  limit,  as  opposed 


o 


to  a  maximum,  of  the  vibrational  frequency  when  the  OH  bond  in  the  linear 


0-H***0  linkage  of  ice  VII  has  its  smallest  length.  From  Fig.  7  this  would 

o 

correspond  to  0.94  A,  which  in  view  of  the  scatter  of  the  data  is  not  much 

°25 

different  from  the  gas  phase  value  of  0.96  A  . 
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X-ray  lattice  parameter,  a,  in  A,  versus  pressure  in  GPa  for  ice  VII  at 
room  temperature. 

Raman  OH-stretching  peak  frequencies  Av,  in  cm  versus  pressure  in  GPa 
for  ice  VII  at  room  temperature. 

Raman  spectra  in  the  OH-stretching  region  for  ice  VII  at  3.2  and  15.5  GPa. 
Raman  versus  x-ray  correlations,  A\>g^  versus  rgg,  for  ice  VII.  Points, 
computed  from  least  squares  equations  for  P  as  a  function  of  a  and  P  as 
a  function  of  AVg^,  see  text.  Curve  A  refers  to  a  linear  and  curve  B 
to  a  cubic  least  squares  fit  to  the  computed  points. 

Mode  Grtineisen  parameter,  Yq^|»  computed  from  method  A  (linear)  and 
method  B  (cubic)  least  squares  fits  of  A\>  versus  rgg,  see  Fig.  4 
and  text.  Points  shown  delineate  region  of  x-ray  and  Raman  data. 

Dashed  lines  refer  to  extrapolations  according  to  method  A  or  method  B. 
Weakly  coupled  oscillator  model  for  the  hydrogen  bonds  in  ice  VII.  The 
0-H  force  constant  is  k^,  and  k2  is  the  force  constant  for  H****0.  The 
figure  defines  relevant  distances  with  d  being  the  proton  displacement 
from  its  equilibrium  position. 

e 

Plot  of  rg^  versus  rgg  both  in  A  for  hydrogen-bonded  systems.  Points 
from  Speakman,  ref.  25.  Solid  line  computed  from  x-ray  and  Raman  data 
for  ice  VII. 
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TABLE  CAPTIONS 

I.  X-ray  lattice  parameter  data  for  ice  VII  to  36  GPa. 

II.  Raman  OH- stretching  peak  frequency  data  for  ice  VII  to  30  GPa. 


o 


o 


Table  I 


p.GPa 

a, A 

2.30 

3.368 

2.30 

3.365 

2.85 

3.337 

4.76 

3.268 

5.45 

3.256 

6.20 

3.231 

8.85 

3.166 

10.07 

3.154 

11.64 

3.109 

11.65 

3.106 

12.47 

3.110 

14.04 

3.078 

15.05 

3.074 

15.58 

3.048 

15.58 

3.045 

17.70 

3.042 

17.95 

3.035 

19.57 

3.029 

20.13 

3.008 

20.32 

3.013 

23.40 

3.011 

24.30 

2.961 

24.57 

2.966 

24.58 

2.968 

26.50 

2.974 

26.50 

2.970 

29.65 

2.940 

30.00 

2.923 

36.00 

2.896 

o 


o 


Table  II 


P.GPa 

-  -1 

Au.Cffl 

p.GPa 

Av.cnf 

2.04 

3308. 

8.10 

3124. 

2.06 

3310. 

8.24 

3102. 

2.06 

3304. 

8.30 

3122. 

2.08 

3306. 

8.57 

3110. 

2.12 

3298. 

8.63 

3102. 

2.18 

3305. 

8.76 

3100. 

2.18 

3310. 

9.50 

3074. 

2.18 

3299. 

9.60 

3082. 

2.26 

3308. 

9.60 

3092. 

2.38 

3308. 

9.90 

3075. 

2.39 

3305. 

10.70 

3050. 

2.49 

3286. 

10.83 

3045. 

2.56 

3288. 

10.89 

3040. 

2.61 

3301. 

10.90 

3044. 

2.62 

3291. 

10.90 

3050. 

2.62 

3288. 

11.00 

3062. 

2.69 

3288. 

11.49 

3028. 

3.02 

3273. 

12.50 

3025. 

3.19 

3263. 

12.67 

3027. 

4.32 

3235. 

13.41 

3009. 

4.44 

3230. 

13.68 

3013. 

4.78 

3216. 

13.96 

3003. 

5.18 

3213. 

14.01 

3001. 

5.68 

3185. 

14.21 

2997. 

5.84 

3180. 

14.38 

2990. 

6.10 

3170. 

14.63 

2994. 

6.37 

3163. 

14.75 

2974. 

6.38 

3155. 

14.99 

2982. 

6.50 

3166. 

15.29 

2975. 

6.80 

3160. 

15.50 

2982. 

7.26 

3147. 

19.30 

2931. 

7.60 

3150. 

22.80 

2893. 

7.70 

3143. 

26.40 

2848. 

7.82 

3109. 

30.00 

2831. 

ICE  VII 
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